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ON THE ROLE OF CYLINDRICAL FUNCTIONS IN 
KANTOROVICH DUALITY 

MARTIJN ZAAL 


Abstract. We study the dual formulation of the Monge-Kantorovich optimal 
transportation problem, in particular under what circumstances it is permitted 
in an infinite dimensional setting to use cylindrical functions, i.e. functions of 
the form (poP where P is a finite-rank operator and is a smooth, compactly 
supported function. In the last section, some examples of applications are 
presented. 


1. Introduction 

Even though optimal transportation problem, first studied by Monge in the paper 
Memoire sur la theorie des deblais et de remblais in 1781 [^, is easy to state, it 
gives rise to rich theory and has proven useful in many applications. Only in 1942, 
Kantorovich provided an answer to the existence question by studying a relaxed 
version of the problem [^. In present time, the Monge-Kantorovich problem is 
understood to be the search for a probability measure 7 on the product X xY with 
prescribed marginals /i and v on X and Y respectively, minimizing the cost 

( 1 . 1 ) / c{x,y)d'y{x,y). 

JxxY 

Since the late 1980’s, the Monge-Kantorovich optimal transport problem has re¬ 
ceived renewed attention and proved to be useful in many fields of mathematics 
including dynamical systems, differential geometry and probability theory. 

For a survey of the Monge-Kantorovich problem and its applications, the reader 
is referred to, for example, m, i or [To] . 

Also introduced by Kantorovich is a problem dnal to the minimization of (HU: 
the maximization of 

(1.2) / (pdy,+ / ipdv 

JX JY 

among all pairs ( 0 , ip) of integrable functions such that (p{x) +ip{y) < c{x, y) every¬ 
where. Being closely related to the Monge-Kantorovich optimal transport problem, 
the dual problem serves as a valuable tool for analyzing problems and questions 
related to the Monge-Kantorovich optimal transport problem. 

The class of functions over which (HU is to be maximized can be reduced consid¬ 
erably (see also O Section 1.1]): the supremum does not change when (p and ip are 
required to be continuous in addition. In case the spaces X and Y are Polish, that 
is, complete separable metric spaces, (p and ip may also be assumed to be bounded. 
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Finally, if X and Y are locally compact, (p and ip may be required to decay to zero. 
One should note, however, that a maximizer of (USD, if it exists, may not possess 
this regularity. 

We will focus on the case where X and Y are vector spaces. The case X = Y = 
E," with a cost function of the form h{x — y) seems to be well-studied (see, for 
instance HI)- In case X or 1" is infinite dimensional, much less seems to be known. 
The Monge-Kantorovich problem when X = is an abstract Wiener with a specific 
cost has been studied by Feyel and Ustiinel [3], but this reference does not explore 
the structure of the dual problem. To the author’s knowledge, no results but the 
general duality results discussed above are known for infinite dimensional vector 
spaces. Here, we will show that, under certain circumstances, the class of functions 
used for the Kantorovich dual problem can be limited to cylindrical functions, more 
precisely functions of the form tp o P, where P is a finite rank operator and (/J is a 
smooth function with compact support. 

The organization of this paper is as follows: in Section [2l some preliminary 
results concerning duality in Polish spaces and Euclidean space are recalled. The 
main results are stated and proven in Section|31 Finally, in Section|31 some examples 
are presented. 


2. Preliminaries 

2.1. General duality. In this section, we will briefly consider the general Kan¬ 
torovich duality theorem. Here and in what follows, we assume that X and Y are 
Polish spaces, and c is a nonnegative, lower semicontinuous function. The value 
- 1-00 is permitted for c. 

Theorem 2.1 (Kantorovich duality). Define the set <I>c to be the set of pairs 
{(j>,ip) € T^(/i) X L^(y) such that 

(2.1) (pix)+ip{y) < c{x,y). 
for y-almost every x and v-almost every y. Then 

( 2 . 2 ) = sup < / (l)dfj,+ / ipdv 

y,Jx Jy 

Moreover, in (12.21) . the supremum may be taken over H [C^{X) x Cl^(Y)). 

See, for instance, O Theorem 1.3]. Actually, going through the proof of this 
theorem, it follows that the supremum can be limited to those (tp, ip) that satisfy 
(EH) everywhere. 

By some carefully chosen modifications of an admissible pair ((p,ip), that is, a 
pair satisfying EH, the supremum in (12.21) can be computed over a smaller class 
of functions. More precisely, a number of additional bounds may be assumed for cp 
and Ip. 

Lemma 2.2. Let $* be the set of all pairs {(p,ip) G C))(X) x C^fY) such that 

(2.3) 

inf (/) = 0 , sup 0 < sup c. Ip > inf {c{x,y) — (p{x)} , swpip>0, 

^ X XxY xeX,yeY Y 

(2.4) inf {c(x, y) - <p(x) - ip(y)} = 0. 

x^X,yGY 
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Then 

(2.5) 


sup < (/)d/.t+ / '0du 

\^Jx Jy 


Proof. Let {4>Y) € C'°(X) x C^{Y) satisfying (12.11) be given, and write 

( 2 . 6 ) V'°( 2 ;) := inf (c(x,.) — ^), A := sup (^ — V'°) + inf 

Y X ^ 

Since c is bounded from below and (j) and ip are bounded, ip^ is bounded from below 
and A is finite. Let then cp := {cp — A)"*". As A > infx cp, we have infx ^ = 0, and 
clearly inixiip^^ — (p) infx Y- By construction, 

(2.7) (p — X < 'ip'^ — inf ip'^. 

As the right hand side is nonnegative, it follows that (p < Y ~ inf a '0'^, that is 
'0'^ — 0 > infx i/)^. We conclude that infjc(0“^ — 0 ) = infx Y- 
Now define 

( 2 . 8 ) 0 := max "I 0 + inf 0 “^, inf <c{x,y) — p{x)\ > 

I X x£X,y£Y I J J 

Again, note that bounds on c, 0 and ip imply that ip is well-defined. Since infx (0'^ — 
0 ) = infx 0 °, we have 

(2.9) c{x, y) — p{x) — ip{y) — inf 0'^ > ip‘^{x) — p{x) — inf 0° > 0 
for any x € X, y gY, and 

(2.10) ^^inf^^ |c(t,2/) - 0(t) - mf^^ - 0(C)}} = 0. 

It follows that 

( 2 . 11 ) inf <1 c(a;,y) - 0 (a;) - 0 (y) !• = 0 . 


x£X 


inf I c{x, y) - p{x) 

X,yGY K 


Given e > 0, there are x G X such that ip‘^{x) < inf a ip‘^ — ^. By (12.71) . (p{x) < |. 
For any y GY such that 


c{x, y) - ip{y) < ip‘'{x) + - <u^ip^ + £ 

A X 


'tpiy) > tp{y) + inf Ip‘' > c{x, y)-£> -£, 

X 


( 2 . 12 ) 

it now follows that 

(2.13) 

that is, supy 0 > 0. Consequently, 

(2.14) 

Finally, 

(2.15) / pdy,+ ipdv> (pdy — X+ ip du + ini ip'^ > / 0d/i- 


< inf {c(.,j/) - 

y^Y I 


< sup c(., y) — sup 'll) < sup c. 

y£Y Y XxY 


fX 


fY 


X 


'll; diy 


fX 


fY 


Since 


(2.16) inf — A = — sup((^ — 'tpc) = inf {c{x, y) — (/)(x) — V’(y)} ^ 0- 

X X x^X,y^Y 

We conclude by applying Theorem 12.11 


□ 
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Remark. Combining the bounds in (12.31) and (12.4p . some more bounds can be de¬ 
rived: 

(2.17) 

— sup c < — supi)) < < sup c, ij) < inf {c{x,.) — ipix)} < inf {c(x ,.)} 

ATxY A XxY a:eX 

whenever (</>,-0) € $*. 

In what follows, it will prove useful to consider a slightly generalized version of 
$*. For 5 > 0, we will denote by the set of pairs (</), -0) S <l>c ("I (C°(X) x C°{Y)) 
such that 

(2.18) —5 < inf (/) < 0, sup(/)< sup c, tp > — snpcj) — 6, suptp>—6. 

^ X XxY X Y 

From this and relaxed versions of (I2.17|) can be derived. 


2.2. Euclidean space. In this section, we will prove that the supremum in the 
Kantorovich dual problem may be restricted to compactly supported smooth func¬ 
tions when X and Y are Euclidean spaces. The proof of the result below is adapted 
from [3 Section 2]. 

As will be come clear in the argument below, some coercivity for the cost function 
is needed. Clearly, it is not realistic to assume that sublevels of the cost c are 
compact: even in case X = K = IR and c{x, y) = \x — y\P for some p > 0, the cost 
is equal to zero on the unbounded set {(a;,^) : x € R}. However, such a cost does 
satisfy the following assumption. 


Assumption 2.3 (Local coercivity). Given any i? > 0, 

0 < (5 < inf sup{c(a:, j/)} 
\x\<Ry^Y 

and M > 0, there exists an i?' > 0 such that |a;| < i? and 


implies \y\ < R'. 


c{x, y) < min 


sup{c(x, T])} — 6,M 

riGY 


Lemma 2.4. Suppose c{x,y) : K,™ x R" — >■ [0, +oo) is lower semi-continuous and 
satisfies Assum,vtion \2.A Then 

(2.19) '^c(/i, i^)= sup ( [ (j){x)dfi+ f 'ip{y)di' 

(0,j/>)e$fn(C“(R”‘)xC“(R")) \Jr^ 

for all fi G tf^{X), V G J^(F) and 6 > 0. 


Proof. Throughout the proof, let y : [0, -l-oo) —)> [0,1] be a smooth function such 
that r]{z) = 1 for z < i and r]{z) = 0 for z > 1, and denote by Br C R* the ball of 
radius r > 0 centered at 0. For fc G IN, p > 0, define 77 ^ G C))“(R^) by 

vM r 

(2.20) Vpix)= afe := / yi\x\)dx. 

akp'^ jRk 

Let then 0 < £ < <5, and suppose that G $*. There exists > 0 be such 

that 


(j>dp < 


e 

6' 


( 2 . 21 ) 
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As (j) is continuous, (j) uniformly on bounded sets as p 4- 0. In particular, 

there exists G (0,1) such that 

( 2 . 22 ) sup ( 77 ™ 

^Rx ” 

Define then 

(2-23) {(t^XBux-Px - §XBr.+pJ , 

where xe is the characteristic function of A, defined by xe{x) = 1 on E and 
Xe(x) = 0 outside E. By construction (j) is smooth, and supported in Bfix+ 2 px- 
Moreover, 

(2.24) 

in Bpix- As (/) < 0 outside B^^, 4> ^ <P everywhere. Finally, everywhere. 

Since sup'!/> > 0, 


(2.25) 


(/)(x) < sup {c(x,y)} . 
yen'* 


By Assumption 12.31 there exists Ry > 0 such that |a;| < Rx and \y\ > Ry implies 
(2.26) 

c{x, y) > min < sup{c(x, y)}, sup ^ > — | > min < ^(x), sup </> > — | > (j){x) — |. 
UeY X } I A J 

Without loss of generality, it may also be assumed that 


(2.27) 


As before, there exists Py G (0,1) such that 


-ipdi/ < 

4 


(2.28) 

Again, defining 
(2.29) 


sup (yZ * ^ - V') 

Brv 


+ e 
< 


^ ■■= vZ * 


X.B Ry - Py 


IXBr 


y+Py 


results in 'i/i < 7 /; on B^y , and i/i < 0 outside B^y. Moreover, i/i > i/’ — | is smooth 
and compactly supported. By construction and (12.261) . 


(2.30) 


^{x) + Ziy) < 


1 

(<j){x)+y{y) - |, 

if 

|x| 

Rxi 

|y| 

< 

Ry, 

< 

1 0(a:) - |> 

if 

|x| 


\y\ 

> 

Ry, 

- 

1 V'(y), 

if 

|x| 

>Rx. 

\y\ 

< 

Ry, 

1 

[o, 

otherwise, 




VI 

i.x,y), 








and (^f), 7 /') G 'hf. Moreover, 


ipdp, ■ 


(2.31) 


t/j du > 


> 


'Br, 


(t>dp - 3 + 


(pdp - 


jBRy 

ipdv — £. 


7/)du — - 


/R" 


An application of Lemma 12.21 now concludes the proof. 


□ 




6 


MARTIJN ZAAL 


3. Frechet space 

3.1. Prechet Spaces and Finite Dimensional Approximations. In this sec¬ 
tion, we focus on the case where X and Y are Frechet spaces. Remember that 
a locally convex topological vector space is called a Frechet space if its topology 
is generated by a translation invariant metric that makes the space complete. In 
the setting of Frechet spaces, one must be careful with the word ‘bounded’: as 
many different metrics may generate the topology, boundedness is not well-defined 
in the usual metric sense. Instead, a subset IF of A is bounded of for every open 
neighborhood C/ of 0, IF C for some t > 0. 

In order to properly define the class of cylindrical functions on a Frechet space A, 
it is necessary to define some kind of coordinate system. One obvious choice would 
be to demand that A has a Schauder basis. It is possible, however, to consider 
a more general class of spaces: we will consider spaces that have the so-called 
approximation property, that is, spaces that admit a family {Pk : A —Aj^gM of 
continuous finite rank operators such that converges to the identity uniformly 
on precompact subsets of A. 

Remark. In fact, it is sufficient to assume that convergence of Pk to the identity 
is pointwise, since pointwise convergence of linear maps from a Frechet space to 
any topological vector space implies uniform convergence on compact subsets. This 
fact is easily proven from the Banach-Steinhaus theorem (Sj 2.5]. It also show that 
{^fe}feeiN is equicontinuous. 

Remark. Note that any Frechet space with the approximation property is also 
separable and hence a Polish space. 

By definition, the range of a finite rank operator is isomorphic to some Euclidean 
space, that is, there exists n S IN and a linear bijection Lp : P\X] —>■ R" that is 
continuous and has continuous inverse. In what follows, we will not make the 
distinction between P and ip o P[X] if there is not chance of confusion. We write 

(3.1) Cyl(A;P):={¥.oP:^eCr(R”)}, 

and call the elements of Cyl(A;P) cylindrical functions with respect to P. A 
function is called cylindrical if it is in Cyl(A; P) for some P. 

In the remainder of this section, we will assume that A and Y are Frechet 
spaces, that {P/dfceK and {Q/cjfceK are families of continuous finite rank operators 
converging to the identity map on A and Y, respectively. 

3.2. Statement and Proof of the Main Results. In order to apply the results 
from Section 12.21 some requirement similar to Assumption 12.31 is needed. More 
precisely, Lemma 12.41 will be applied to the cost composed with Pk and Qk, which 
means that this function has to satisfy Assumption [531 One can make the following 
assumption on c itself to ensure that this is the case. 

Assumption 3.1. Given any bounded V C A, 0 < d < inf 2 ,gB^( 3 ;) supygY{c(a^: J/)} 
and M > 0, there exists a bounded set IF C F such that x gV and 

c{x, y) < min < sup{c(a;, y)} — 6, M 
ir&Y 


implies y S IF. 
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In what follows, we will not assume that this assumption holds, but use a logically 
weaker requirement. We will see later, though, that Assumption 13.11 holds for some 
interesting examples. 

Theorem 3.2. For fc G IN, define Ck '■ X xY —>■ [0, +oo] by Ck{x, y) = c{PkX, QkV), 
and suppose that the restriction of c to Pfe[A] x (5fe[T] satisfies Assumption IS.M 
Then 

(3.2) < liminf sup < / (/)d/i+ / i/>dul 

(0,i/')e<I>f,^n(Cyl(A;Pfe)xCyl(Y;Qfc) UA Jv J 

for any y G J3^(X), v G IX^fY) and any i5 > 0. 

Proof. Let e > 0, and choose ((jijif) G such that 

(3.3) f 4><i.y+ [ tpdi' + S. 

Jx Jy 

By assumption, (p o Pk ^ (f, and i/> o Qk if. As p and if are bounded, we obtain 

(3.4) ‘^c(M: < lim if 4> o Pk dfi + [ -if oQk du ) + e. 

fe-^-oo \JX JY / 

Since {(j), ip) satisfies (|2.ip for all a: G A, y G Y, the restrictions of (p and ip to Pk[X] 
and Qk\Y\, respectively, satisfy (12.ip as well. Therefore, Lemma [2.41 implies 


(3.5) / (poPkdp.+ / ipoQkdv 

Jx Jy 

= / 0d(Pfe)#^+ / ipd{Qk)i).v 

jR'” jR" 


< sup 

(Ab)e'i>^n(C~(R™)xC~(R")) tJR 


I (^d(Pfe)#/r+ / ipd(Qk)n.i 
Jr™ Jy 


Now (p := (p o Pk G Cyl(A;Pfe) and ip := ip o Qk & Cyl(T;Qfc) whenever (p G 
C)N(R'"), Ip G C)N(R”). Moreover, (p,ip) G 4)^^ by construction. Therefore 


(3.6) / po Pkdfj,+ / ipoQkdu 

Jx Jy 

< sup 

(0,^)e$«^n(Cyl(X;Pfc)xCyl(Y:Qfc)) 

Combining this with (13.4p results in 




(3.7) '^ciPii') Fidmird sup 

(<^,Vi)e<I-^^n(Cyl(A;Pfc)xCyl(y:Qfc)) 



pdy + 



+ e. 


As e > 0 is arbitrary, we conclude. 


□ 


In order to obtain a lower bound, more structure has to be assumed for c. One 
possibility is to assume that c is continuous and its growth can be controlled in 
a certain way. Under these assumptions, the converse of (13.211 can be shown by 
considering the Monge-Kantorovich optimal transport problem, rather than the 
dual problem. 
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Proposition 3.3. Suppose that, additionally, c is upper semicontinuous (and hence 
continuous), and there exists nondecreasing functions /, g : [0, +oo) —>■ [0, +oo) such 
that 

(3.8) c{Xx,y) < f{X)c{x,y), c{x,Xy) < g{X)c{x,y), 

for any X > 0, x € X and y GY. Then 

(3.9) 

‘^c(/J., J^)'^c(M) ^) ^ l™sup sup < / (j)dp,+ / tpdi'? 

fe->.oo (,^,V<)e$^j^n(Cyl(A;Pfc)xCyl(Y;Qfc) Ux Jy J 

Proof. In case ‘^c(Mi '^) = +oo there is nothing to prove, so assume that 7 G r(/i, v) 
is such that 

(3.10) / c{x,y)d'^{x,y) <+ 00 . 

JxxY 

By assumption, c{PkX,Qky) < f{X)g{X)c{x,y), which means that 


c{x,y)d-f{x,y) > 


IXxY 


/ limsupc(Pfea:, ( 5 fcy)d 7 (x, y) 
'xxY k^oo 


(3.11) 


>limsup/ c{PkX,Qky)d'y{x,y) 

k—^co JX xY 

> limsup'^c ((P/c)#y, {Qk)n.v) 

k—¥oo 


As Pk [A'] and Qk \Y] isometrically isomorphic to some finite-dimensional Euclidean 
space, the right hand side is equal to the right hand side of (13.9p by Theorem 12.41 
Taking the infimum over 7 G r(/i, v) now implies the result. □ 


Combining the previous two results implies that the optimal transport cost can 
be written as the limit of suprema over classes of cylindrical functions. However, 
if the operators Pk and Qk can be chosen such that Ck is in fact at most c, it is no 
longer needed to take a limit of suprema. 


Theorem 3.4. In case c{PkX, Qkv) < c[x, y) for allk gM, x G X and y GY, 

(3.12) = sup < / (j)dpL-G / ipdix 

fceK,(0,V')e$fn(Cyl(X;Pfc)xCyl(Y;Qfc) Ux Jy 

whenever p G tX^(X), v G I^iY) and ^ > 0. 

Proof. By assumption, C all fc G IN, 5 > 0. Therefore, 

'^c(/i, p)= sup < / (j)dp+ ipdix 
Ux Jy 


(3.13) 


> sup 

fceK,(0,V')e$fn(Cyl(A;Pfc)xCyl(Y;(3fc) Ux 


cj)dp + Tpdv 


> sup 


/ceK.(0,V-)e$fj^n(Cyl(X;Pfc)xCyl(Y;Q*,) Ux 

By Theorem 13.21 this implies the result. 


(j)dp+ ifdv 


lY 


□ 


Remark. Note that the proof of Proposition 13.31 also works under the assumptions 
of this theorem. Therefore, (1^ also holds in this case. 
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4. Some examples 

4.1. Separable Hilbert Space with Translation Invariant Cost. Let us con¬ 
sider the case X = Y = H some separable Hilbert space with inner product (.,.) 
and norm |.|, and c{x,y) = h{\x — y\) for some lower semicontinuous function 
h : [ 0 , -boo) —>■ [ 0 , -boo] such that h{z) sup[Q__|_j,o) h as z ^ -boo. 

It is not very difficult to check that Assumption 13.11 holds in this case. Given 
M < sup[Q_^oo) h, the cost c{x,y) can only be less than M ii x — y is restricted to 
some bounded set. If x is already restricted to a bounded set, this means that y 
must also be restricted to a bounded set. Therefore, Theorem o can be applied 
in this case. 

The maps Pk and Qk can be constructed from a complete orthonormal sequence 

{ejljgiN: 

k 

(4.1) PkX = QkX = '^{x, ej)ej 

1=1 

If, additionally, h is continuous. Proposition 13.31 can be applied as well. There¬ 
fore, we conclude that 

(4.2) '^c(/i, u)= lim sup I (j)dfj,+ i/jdub. 

(<b,b')e4'?^^n(Cyl(A;Pfc)xCyl(Y;Qfe) [Jx Jy J 

If we assume h to be nondecreasing instead, the example fits the requirements 
of Theorem 13.41 and we obtain (13.121) . 

It is easy to see that the similar results hold if 

(4.3) c(x,y) = h 

with hj > 0 lower semicontinuous with hj{z) —1 -boo as z —1 -boo for every j € IN, 
and h continuous and increasing. Although c may not satisfy Assumption 13.11 it is 
obvious that 

(4.4) ck{x,y) = h ^^/ij((a;,ej))j 

satisfies 12.31 for every k. Depending on the continuity, growth and summability of 
the hj, Proposition 13.31 may or may not be applicable in this case. If, however, all 
hj are increasing, Theorem 13.41 can be applied. 

4.2. Banach Spaces With a Schauder Base. If A = T is some Banach space 
that has a Schauder base, there is a canonical map Pk- However, since it is not a 
priori true that ||Pfc|| < 1, Theorem 13.41 cannot be applied a priori, even if c(x,y) = 
jja; — yp for some p > 0 . 

If c has a structure similar to the structure in the previous section, it is possible 
to apply Theorem 13.21 and Proposition 13.31 to obtain (14.21) . 

4.3. Wiener Space. An important example of an infinite dimensional vector space 
where optimal transportation may be of interest is Wiener space (see also [3] . If IT 
is a separable Frechet space, and 7 is a Gaussian measure with support IT, then 
the Cameron-Martin space iJ of 7 is dense in IT and separable (see, for example. 
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m- Therefore, there exists a complete orthonormal sequence {hfc}fcg]N in H such 
that the functionals {hk, ■)h extend by continuity to W. For such a sequence, the 
orthogonal projections Pk onto the span Fk of {hi,..., hk} defined by 


k 



(4.5) 


also extend by continuity to W. Moreover, as the sequence {hk} is complete and 
H is dense in W, the union of Fk is dense in FI. However, it is not in general true 
that PkX ^ X'm all of IT as {Pfe} is not necessarily equicontinuous as a collection 
of maps from W to W. 

In case Pk does converge pointwise to the identity, the results from the previous 
section can be applied. Particularly interesting is the cost function that |a; — y\H 
whenever x — y G H and +oo otherwise. It is not difficult to show that this cost 
function is lower semicontinuous and that c{PkX,Pky) < c{x,y) for all x,y G W. 
Therefore, Theorem 13.41 can be applied in this case. 
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